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ABSTRACT 

Values of black hole masses are frequently determined with the help of the rever- 
beration method. This method requires a specific geometrical factor related to the 
distribution of the orbits of the Broad Line Region clouds. Onken et al. determined 
the value f 2 = 1.37 ± 0.45 from the black hole mass - dispersion relation. In this 
paper we determine this factor using an independent mass determination from the 
X-ray variance method for a number of Seyfert 1 galaxies and comparing them with 
the reverberation results by Peterson et al. We obtain mean value f 2 = 1.12 ± 0.54, 
consistent with Onken et al. Both values are larger than the value 0.75 corresponding 
to a spherical geometry. It indicates that most probably all values of the black hole 
masses obtained with the use of the Kaspi et al. formulae should be multiplied by a 
factor of ~ 1.7. This also shows that the Broad Line Region is rather flat, and hints 
for a dependence of the factor f 2 on a source inclination seem to be present in the 
data. 
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1 INTRODUCTION 

The value of the black hole mass is one of the key parameters 
in the description of the accretion process in an active galac- 
tic nucleus (AGN). Much effort was devoted to develop re- 
liable and efficient methods of mass determination from the 
observational data. Current popular methods include the 
reverberation mapping (Blandford & McKee 1982; Wandel, 
Peterson & Malkan 1999; Kaspi et al. 2000; Onken & Peter- 
son 2002; Vestergaard 2002), the stellar and gas kinematics 
(Kormendy & Richstone 1995; Nelson & Whittle 1995; Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000; Woo & Urry 
2002; Tremaine et al 2002; Verolme et al. 2002; Peterson et 
al. 2004; Onken et al. 2004), determination via water maser 
emission (Miyoshi et al. 1995; Greenhill et al. 2003), the 
method based on BH mass - bulge mass or bulge - luminos- 
ity relation (Magorrian et al. 1998; Kormendy & Gebhardt 
2001; McLure & Dunlop 2002; Haring & Rix 2004), based 
on the disk luminosity estimation (Tripp, Bechtold & Green 
1994; Collin et al. 2002) or just the H a line analysis (Green & 
Ho 2005). These methods were recently supplemented by the 
methods based on X-ray variability: power spectrum break 
(Papadakis 2004; McHardy et al. 2005) and excess variance 
(Nikolajuk et al. 2004; O'Neill et al. 2005). 



Each of the methods require the knowledge of a cer- 
tain normalization factor, and is biased by some systematic 
errors, so the comparison of two independent methods is 
extremely important. 

In the present paper we compare the results of the re- 
verberation approach with the results of the X-ray excess 
variance for several Seyfert 1 galaxies. The methods of mass 
determination and the approach to their comparison is given 
in Section [5] the results for mean normalization factors and 
the dependence on the source inclination are shown in Sec- 
tion [3] and we discuss the results in Section 2] 



2 METHOD 

2.1 Selected sample of Seyfert galaxies 

We consider a sample of Seyfert 1 galaxies which were both 
a subject of optical and X-ray monitoring (see Table 0. 
We do not include Narrow Line Seyfert 1 galaxies since 
there are many indications of peculiar properties of these 
sources in comparison to normal Seyfert 1 objects (smaller 
black hole mass to bulge mass ratio, Wandel 1999; Mathur, 
Kuraszkiewicz & Czerny 2001, higher variability amplitude, 
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Bian & Zhao 2003; Nikoiajuk, Papadakis & Czerny 2004; 
Markowitz & Edelson 2004). In particular, we do not in- 
clude NGC 4395 since it might be a NLS1 galaxy (FWHM 
of U/3 line 1500 km s -1 ; Kraemer et al. 1999). This leaves 
us with 13 sources with black hole masses in the range 
~ 10 7 - 3 x 1O 8 M . 

We use the reverberation results of Peterson et al. 
(2004) since they are based on the most advanced analy- 
sis of the line profiles. Most of the X-ray variability data 
were taken from RXTE and ASCA databases (see Nikoia- 
juk et al. 2004). For F9 we take the excess variances and the 
durations of all 7 observations from Turner et al. (1999). 
Turner et al. give both the full duration of the observation 
and the effective time actually covered by the data. We use 
the full duration To as T since gaps in the data affect less 
the variance than the extension of the monitoring. This is 
related to the steep power spectrum and the dominance of 
the lower frequencies in variability. 

Measurements of the inclination angles of the sources 
were taken mostly from Nandra et al. (1997b). For IC 
4329A, NGC 5548 and NGC 4593 we took inclinations from 
Mushotzky et al. (1995) and Guainazzi et al. (1999), corre- 
spondingly. These measurements are based on the shape of 
the iron Fe Ka line. For 3C 120 and 3C 390.3 we adopted 
the values estimated by Eracleous & Halpern (1998) and 
Ballantyne et al. (2004) from observations of radio jets. 

2.2 Black hole mass from the reverberation 
method 

Reverberation method allowed to measure the mass directly 
in several Seyfert galaxies (Wandel et al. 1999; Kaspi et 
al. 2000). The basic formula underlying the reverberation 
method reads 



M = f 



2 "FWHM^BLR 



G 



(1) 



where M is the black hole mass, G - the gravitational con- 
stant, -Rblr - the radius of the Broad Line Region, wfwhm - 
the Full Width Half Maximum (FWHM) of H/3 line and f 
is the squared geometrical factor describing the distribution 
of the orbits of Broad Line Region (hereafter BLR) clouds. 
In a recent more sophisticated approach of Peterson et al. 
(2004) to the line profile analysis the «fwem was replaced 
with the dispersion o"ii ne : 



M = U 



® line CTccnt 



(2) 



where r ccnt is the centroid of the cross-correlation function, 
c is the light velocity. In the sample of Peterson et al. sta- 
tistically there is a relation between the two factors 

f « \U ■ (3) 

Onken et al. (2004) determined the scaling coefficient 
in equation @ from the correlation between the black 
hole mass and the bulge/spheroid stellar velocity disper- 
sion. Their value of / CT = 5.5 ± 1.8 was adopted by Peter- 
son et al. (2004). We use their results in our inclination- 
independent approach. In this approach there is no space 
for any anisotropy of the BLR. 

However, BLR is unlikely to be spherically symmetric 
(e.g. Done & Krolik 1996; Krolik 2001, Collin & Kawaguchi 



2004). Most plausible geometry, particularly for the Low 
Ionization Line Region, is a disk-like wind (Collin-Souffrin et 
al. 1986; Chiang & Murray 1996; Hutchings et al. 2001; Kol- 
latschny 2003) . Unfortunately, the theoretical predictions of 
the shape of an emission line is quite sensitive to the as- 
sumptions about the optical depth of the wind and the dis- 
tribution of the emissivity (Murray & Chiang 1997). 

Therefore, we consider, as a possibility, that the BLR 
velocity may be represented by a combination of a random 
isotropic component, v T , and a component only in plane of 
the disc, v p (Wills & Browne 1986; McLure & Dunlop 2001). 
Therefore, the masses of the black hole in AGNs can be 
represented by 



M = 



"fwhm-Rblh 



4(£ 2 + sin 2 i) 



G 



(4) 



where £ = v r /v p . It differs by a factor of two in front of sin i 
from the formula of Krolik (2001) and of Collin & Kawaguchi 
(2004). 

This relation means that / CT is not an inclination- 
independent coefficient but 



^ + sin" 



(5) 



If f = 0, the value / CT = 5.5 determined by Onken et al. (2004) 
corresponds to a representative inclination angle i = 25.2°. 

The use of this formula requires the knowledge of £ and i 
for every object separately, and the exact dependence on the 
inclination angle is not well justified. However, it is useful 
for search of any possible traces of anisotropy. 



2.3 Black hole mass from the X-ray variability 

We determine the black hole mass using the excess variance 
method. The method is based on assumptions that (i) the 
high frequency tail of the power spectrum has a slope of -2 
(ii) the high frequency break scales inversely with the black 
hole mass (iii) the value of the power times frequency at 
the high frequency break is universal for all objects, inde- 
pendent from mass. In this paper we improve slightly this 
method, in comparison to the method described by Nikoia- 
juk et al. (2004). We change: (i) the method of obtaining 
the final value of A/bh for a given object from results for 
individual lightcurves (ii) the method of calculation of the 
errors (iii) the value of the normalization constant present 
in our method. 

The X-ray variance method (Nikoiajuk et al. 2004) is 
based on the scaling of the normalized variance of X-ray 
light curves, o"n XS , with black hole mass 

T — 2At 

(T 

nxs 

Here T is the duration of a single X-ray light curve in seconds 
and At is its bin size in seconds, as well. a^ a is in units of 
(rms/mean) 2 . T cannot be longer than the time scale defined 
by the high frequency break of the power spectrum. Longer 
light curves can be chopped and used to determine several 
independent values of the variance. The normalized excess 
variance, o"n XS , is defined as in Nandra et al. (1997a) and 
Turner et al. (1999). 

From each light curve of a given object we obtain the 
fcth individual black hole masses A/bh,*, using equation ©. 
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Table 1. The sources used in our study. 



Name 


Type 


n,f Peters 
M BH,<t2 


M BH, CT 2 


r = 


7i /Peters 


log L x 


r 3 


-10 




^Nandra 


^oth. papers (^af.) 


i 








(10 7 M Q ) 


(10 7 Mq) 














(eV) 


(dag) 


(deg) 


(deg) 


3C 120 


SI 


.5 


cc +3 .14 
5.55_ 2 2g 


- „ —-1-90 d 

±y " ' -8.9 


3 


.55 


± 3.11 


43.95t 


1 


.79 


70 


60+1° 


< 14( fc )(l,2) 




3C 390.3 


SI 


.5 


28.7 ±6.4 


o-f-14 8 

22 -3I 8 4 5 8 





.78 


± 0.44 


44. 20* 










19 < i < 33( b '(2) 


_ „_i_i n 

29lg 


Ark 120 


SI 


.0 


15.0 ± 1.9 


lO.St^ 3 





.72 


± 0.72 


43.88 












-. .-.4-9(1 

30Ii" 


IC 4329A 


C 1 

SI 


.2 


< 2.78 


r ,ir»+3.62 

5 - 42 Il.9g 




> 


1.95 


43.59 


1 


.70 


80 


_i_15 
22 -22 


< 24' a J (3) 


< 18 


Mrk 509 


SI 


.5 


14.3 ± 1.2 


8.00+*;™ 





.56 


±0.23 


44.03 


1 


.76 


70 


a 1 +49 




35+ 9 


NGC 3227 


SI 


.5 


A. 22 ± 2.14 


, i i +1.76 
4 - n -1.17 





.97 


±0.60 


41.66 


1 


.52 


120 


20 -10 




nc+9 
26 g 


NGC 3516 


SI 


.5 


4.27 ± 1.46 


2 95+ 160 





.69 


±0.37 


43.08 


1 


.74 


90 


271* 




31+]° 


NGC 3783 


SI 


.5 


2.98 ±0.54 


r, iq+0.92 

z - lo -0.70 





.71 


±0.30 


42.90 


1 


.54 


115 


26+ 9 
ZD -10 




30+ 7 


NGC 4151 


SI 


.5 


1.33 ±0.46 


2 92+ - 97 


2 


.19 


±0.98 


42.62 


1 


.55 


100 


nq + 12 
ZO -13 


32 + 23(a) (4) 




NGC 4593 


SI 


.0 


< 1.47 


67+ - 29 
u -°'-0.45 




> 


0.45 


42.98 


1 


.78 


90 


45 +45 
40 -45 


< 39 


NGC 5548 


SI 


.5 


6.71 ±0.26 


14 1+ 211 
4.1 


2 


.10 


± 1.88 


43.41 


1 


.75 


90 


47+43 

*'-41 


15 < i < 38< a )(3) 


1711 


NGC 7469 


SI 


.5 


1.22 ±0.14 


9 co+0.65 
2 - 58 -0.51 


2 


.11 


±0.53 


43.25 


1 


.78 


130 


19+71 
iy -19 




17± 2 2 


F 9 


SI 


.2 


25.5 ±5.6 


8 44+ 4 - 78 





.33 


±0.15 


43.91 


1 


.83 


120 


4 5 +45 




48±i| 



Col. (1) lists the object name and type (S1.0-S1.5 denote Seyferts galaxies). Col. (3) shows the masses of black holes obtained from 
reverberation method by Peterson et al. (2004), who take the geometrical factor {/ CT } = 5.5. Col. (4) shows the black holes masses 
obtained from the variance method. The ratios of the two masses are given in Col. (5). Col. (6) Log of 2-10 keV luminosity in units of 
erg s — 1 . All luminosities are given for: Ho = 75 km s — 1 and go = 0.5. The values with the superscript * were taken from Green et al. 

(1993) and transformed to this cosmology. The other values were directly taken from O'Neill et al. (2005). In Col. (7) and (8) are 
power-law index and equivalent width of Fe Kq line, respectively. All values were taken from Nandra et al. (1997b). Col. (9) shows the 
inclination angles given by Nandra et al. (1997b). Col. (10) shows the same quantity but taken from other papers. The superscripts ( a > 
and ( 6 ) in Col. (10) denote inclination obtained from observation of the Fe Ka line and the radio jet correspondingly. The numbers in 
parentheses correspond to the following references: (1) Ballantyne et al. (2004), (2) Eracleous & Halpern (1998), (3) Mushotzky et al. 
(1995), (4) Guainazzi et al. (1999). Col. (11) shows the inclination angle calculated from equation 1121 under the assumptions £ = 0. 



In order to calculate the final value of Mbh from several 
lightcurves, we fit this set of individual Mbh,*- For each 
lightcurve we calculate the coefficient At 



(T - 2At) h 



(7) 



and the final coefficient A is calculated as a weighted mean 
from the values of Ak, with the weight coefficient given by 
\(T-2M) k f, i.e. 



Efi c i(^x 5 ) fc (r-2At) fc 



A = 



(8) 



(9) 



Here (T — 2At)^ and (o"n Xa )fe are results from several mea- 
surements. The number of the lightcurves is iVic- 

The final black hole mass in the object is derived di- 
rectly from the obtained value of A and equation © 

»„ c 
M — — . 

A 

The value of the constant C in equation © must be 
determined by applying the method to the source with a 
known mass. Cyg X-l was selected as a reference in the work 
of Nikolajuk et al. (2004). The mass of the black hole in Cyg 
X-l was assumed to be 10 Mq and the value of the constant 
C = 0.96±0.02 Mq s" 1 was derived. However, recent careful 
study of the evolutionary history of this binary star showed 
that the most likely value for the mass of the black hole in 
the system is 20 ± 5 Mq. Therefore, we assume 20 Mq as 
the best estimate of the black hole mass in Cyg X-l, and 
the constant C in equation JHJ is consequently two times 



higher: C = 1.92 ± 0.5 Mq s~\ The error of this value is 
predominantly connected with uncertainty of the Cyg X-l 
black hole mass. 

The accuracy of the black hole mass determination in 
the X-ray variability method depends strongly on the ef- 
fect of statistical error of variance measurement and power 
leaking from long timescales. Therefore, in our analysis we 
estimate the error by performing Monte Carlo simulations 
for each source separately. A few hundred sets of artificial 
data were generated for each of the sources, using the Tim- 
mer & Konig (1995) algorithm, exponent of the lightcurve 
was calculated to account for the log-normal distribution 
(see Uttley, McHardy & Vaughan 2005), and each set was 
analyzed like the data, thus giving a distribution of the val- 
ues of A, so its error was determined from the dispersion of 
the resulting distribution. The mass error was determined 
from the error of A. 

The error in the normalization constant was discussed 
separately since the results of Peterson et al. (2004) also did 
not include the coefficient error. 



2.4 Comparison of the two methods and search 
for inclination effects 

Assuming the absence of any inclination-dependent trend in 
BLR by calculating the ratio r of the two mass measure- 
ments for all objects in our sample separately: 



]\/f Peters 
JW BH,rcv 



(10) 
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where M BH a i is the black hole mass measured by the excess 
variance method and MjJ'™ is taken from Peterson et al. 
(2004). We next calculate the mean and the median values 
from the obtained distribution. If the obtained mean value 
(r) is consistent with 1 within observational error, the two 
methods statistically give consistent results. 

Next we allow for an anisotropy of the BLR. We as- 
sume that the excess variance method is independent on 
the source inclination but for the reverberation results we 
adopt the coefficient f a given by equation © instead of a 
fixed factor of 5.5. Under this assumption the ratio r of the 
two measurements should follow a relation 



5.5(£ 2 + sin 2 i) ' V ; 

We check the possible presence of such trend using an inde- 
pendent measurements of the sources inclinations (see Sec- 
tion Eg. 

We also show the results of a complementary approach, 
following the general method of Wu & Han (2001) . We adopt 
the relation given by equation 1111 . we determine the incli- 
nation angle of each object from the formula 



i = arcsin y' — - £ 2 , (12) 

and we compare them to independent measurements of the 
inclination. Errors of the calculated inclination angle are 
based on the error propagation theory (Bevington & Robin- 
son, 1969). 



3 RESULTS 

3.1 Mass measurements 

The values of the black hole masses, M^Tcvi determined 
by Peterson et al. (2004) with = 5.5 are given in the 
third column of Table 1. In fourth column we give the masses 
M BH CT 2 determined with the excess variance method, as de- 
scribed in Section ^..'-!! Generally, the two measurements cor- 
relate quite nicely, as shown in Fig. and Fig. [5] 

Statistical measurement errors are usually smaller in 
the case of reverberation method. On the other hand, in 
two cases (IC 4329A and NGC 4593) reverberation method 
gave only upper limits while the excess variance method gave 
both the upper and lower limits, consistent within two sigma 
error with the results from the velocity dispersion method. 

In fifth column of Table 1 we give the mass ratio 
r — M BHiCT 2/M B H t j iev- The distribution shows significant 
scatter, with the minimum value obtained for Fairall 9 and 
the maximum value obtained for 3C 120. Therefore, we cal- 
culate the mean and the median value. In calculations of the 
weighted mean we exclude IC 4329A and NGC 4593, because 
of the upper limits in the determination of M B H,rev from 
those galaxies. The weighted mean is equal to 0.81 ± 0.15. 
The given error is only statistical one sigma error. This in- 
dicates that both methods are normalized rather accurately 
and there is no strong need to introduce the systematic er- 
rors due to f a and C constants. This is very encouraging 
since the level of disagreement is much lower than allowed 
for by systematic errors entering the two methods as er- 
rors of the normalization constants: / CT = 5.5(1 ± 0.33) and 




Log M BH from reverberation [M ] 



Figure 1. Comparison of the black hole masses, M BH CT 2, taken 
from the X-ray variance method with the masses, M B h revi taken 
from the reverberation method. The solid line shows the expected 
M BH (T 2 = M BH ,rev relation. 




Log M BH from velocity dispersion [M Q ] 



Figure 2. The black hole masses obtained from the reverbera- 
tion method versus the masses of the same object obtained from 
the M B h-o* relation. The values of the M B n,rev are taken from 
Peterson et al. (2004). The Msj^disp values of the velocity dis- 
persion method are calculated using o> values obtained by Onken 
et al. (2004) and the formula of Tremaine et al. (2002). The solid 
line shows the equality of the values M BB ,rev and M B jj,disp- 
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Figure 3. Mass ratio r versus inclination taken from literature 
('Nandra and *oth. papers! see references in Table 1). Continuous 
line shows the best fit to the data under assumption that £ = 0.19. 
Dotted line represents the same fit but for £ = 0.0 (see text for 
details). 



Spearman rank correlation test and in both cases we do not 
confirm any correlation. In the case of r-r3_io k c v a Spear- 
man Correlation Coefficient R = —0.06. In the second case, 
where we test t-Wkq once again the coefficient R does not 
show any correlation (R = —0.07). 

Uttley and McHardy (2005), McHardy et al. (2005) 
have suggested that the deviation from linear scaling be- 
tween the high frequency break and Mbh may by a function 
of accretion rate. In order to check this hypothesis and its 
possible influence to our mass estimation, A/ BH ct 2, we ex- 
amine a new correlation. We check whether r correlate with 
the ratio Lx /M BH a i (i.e. the ratio of 2-10 keV luminosities 
divided by M BH a i~). We use sample of eleven objects (all 
those objects, which have calculated r) and find an appro- 
priate values of Lx in papers of O'Neill et al. (2005) and 
Green et al. (1993). The cosmological parameters which we 
use are: Ho = 75 km s _1 and go = 0.5. The Spearman Cor- 
relation Coefficient R, which we obtained, is equal to -0.54. 
We can say that r probably slightly anticorrelates with the 
ratio tx/M BH (T 2. Nevertheless, it is difficult clearly confirm 
this hypothesis, because only \R\ > 0.73 (for the sample of 
11th pairs) indicates the presence of the correlation at the 
95% level of confidence. 

Assuming a BLR model outlined in Section [2.41 we ex- 
pect the inclinations to follow the equation 11211 . Since £ 2 
in equation 1121 is unknown, we attempted to determine it 
from the requirement of the best fit to the data. We searched 
for a minimum of the function: 



C = 1.9(1 ±0.26) Mq s _1 . It also means that indeed statis- 
tical errors are likely to be dominating individual measure- 
ments in both methods. 

No apparent trend with the black hole mass is present, 
so probably no large systematic error related to the black 
hole mass is involved. 



3.2 Search for BLR anisotropy 

If the BLR is anisotropic this anisotropy is likely to show up 
in a different way in the two black hole mass measurement 
methods. Therefore, in Fig.|H]we show the relation between 
the mass ratio, r, and the inclination angles measured in- 
dependently (see Section [2. in . Errors are large but the two 
quantities seem to be correlated: larger mass measurement 
ratios are more likely to be derived for lower source incli- 
nation. The Spearman rank correlation test for all ten pairs 
does not confirm this impression (R — 0.633). However, if 
the inclination angle of NGC 5548 (47°) from Nandra et al. 
(1997b) is replaced with the value of 26.5° (mean value from 
Mushotzky et al. 1995, see column (10) of Table 1, also from 
the iron line profile), the presence of the correlation is con- 
firmed (R — 0.85). Without this object, the correlation is 
excellent (R = 0.975). 

In order to investigate an importance of our correla- 
tion we check up if the mass ratio correlated with power law 
index, F3_io kcv, or with line equivalent width, Wkc- We 
take the appropriate values in Nandra et al. (1997b), and 
we include them in Table Q If we do not find any corre- 
lation between F3_io kcv and r and/or between Wko and 
r then our correlation, the mass ratio-inclination angles, is 
more believable. For the sample of our objects we do the 



(5i) 2 



(13) 



where (5i) 2 included both the measurement error of the in- 
clination angle (given in Table 1) and the error in prediction 
of ianai from equation 1121 due to an error in the measure- 
ment of r. Resulting \ 2 ls lowest for £ = 0.19 but the errors 
are so large that £ = 0.0 is also acceptable. We show both 
fits to the data in Fig. |3] 

One object (3C 120) shows significant departure from 
the overall mean. High value of the r is likely to be due to an 
extreme value of the viewing angle. Small inclination angle 
i decreases the line widths (see equation^. 3C 120 is most 
probably viewed face on since it shows superluminal motion 
(e.g. Homan et al. 2001). 

In Fig. 2] we compare the determination of the incli- 
nation from equation 112H (Table 1, last column) with the 
results (mostly) from the shape of the iron line (see Sec- 
tion 123- The agreement between the two methods is en- 
couraging. It both supports the existence of the anisotropic 
part in the velocity field of BLR and the correctness of the 
inclination determination by Nandra et al. (1997b). Recent 
results from XMM-Newton telescope frequently do not sup- 
port the existence of the broad relativistically smeared iron 
line (e.g. Mkn 509, Page et al. 2003; NGC 3227, Gondoin 
et al. 2003; NGC 4151, Schurch et al. 2003; NGC 5548, 
Pounds et al. 2003; NGC 3783, Reeves et al. 2004; NGC 
4593, Reynolds et al. 2004; NGC 3516, Turner et al. 2005). 
It is possible that high quality data is sensitive to the de- 
partures of the actual line from the fitted shape and reject 
the presence of such feature. This departure may be, for 
example, related to the warm absorber effect since narrow 
absorption features superimposed onto a broad line mod- 



© 0000 RAS, MNRAS 000, 000-000 



6 Nikoiajuk, Czerny, Ziolkowski, Gierlinski 




20 40 60 

inclination our [cleg] 



Figure 4. Comparison of the inclination angles of the accretion 
discs obtained from the X-ray variance method with the incli- 
nation angles taken from different papers. The inclinations from 
literature are derived mainly from Fe Ka observations. Other val- 
ues are based on the radio jet observations. The references of the 
papers are included in Col. (10) of the Table 1. The solid line 

shows the 2 Irom literature — *from the X — ray variance relation. Dif- 
ferent lines, which indicate errors of i are only to avoid mistakes 
and 'well' viewing. 

ify the shape of the line making it more symmetric, almost 
two-Gaussian (Rozariska et al. 2006). 

We were also able to determine the inclination for Ark 
120 (no previous inclination measurements were done for 
this object). Measured value is quite reasonable for a Seyfert 
1 galaxy. It will be interesting to see whether future inde- 
pendent determinations will confirm our result. 



4 DISCUSSION 

4.1 Consistency of black hole mass measurement 
methods 

We compare two methods of black hole mass determination: 
reverberation method and X-ray excess variance method. 
Both methods include important scaling factors. 

Reverberation method is based on determination of the 
width of the broad lines and the delay in their response 
to variable continuum. It requires an assumption about the 
BLR geometry, in a form of a specific value of the geometri- 
cal factor. Onken et al. (2004) obtained this factor by using 
the black hole mass measurements derived from the velocity 
dispersion measurements. This factor was used by Peterson 
et al. (2004) for the reverberation measurements used in our 
paper. 

Excess variance method is based on a scaling of the 
X-ray power spectra with black hole mass. It requires an 
a priori knowledge of black hole mass in a reference object. 



It is convenient to use the Galactic X-ray source Cyg X-l in 
its hard state to scale the method. The best current deter- 
mination of the black hole mass in this object is 20 ± 5Mq 
(Ziolkowski 2005), and we adopt this value in the present 
analysis. 

We also attempted to use another Galactic source, 
GX 339-4, as a reference. We have extracted power spectra 
from the hard state during the rise of the 2002/2003 out- 
burst (observations 0-24 from Belloni et al. 2005) from the 
PCA onboard RXTE. We have corrected the spectra for the 
Poisson noise and dead time effects of the instrument. Then, 
we integrated them between 10 and 128 Hz in order to deter- 
mine the constant C. We have assumed the mass of 10 Mq 
(Zdziarski et al. 2004) and found C GX 339-4 = (0.91 ± 0.04) 
Mq s~ , by averaging over all observations. It is interest- 
ing to note that the value of C for individual observations 
remained roughly constant throughout the wide range of ob- 
served luminosities (count rate ~400-900 s _1 per one PCA 
detector). The resulting constant C is less than the value of 
C C yg x-i = (1.92 ± 0.5) M s~\ found in Section IO but 
uncertainties due to mass estimation are large. In particular, 
the mass of GX 339-4 is not known very well, with only the 
lower limit from the mass function, 5.8 ± O.5M0 (Hynes et 
al. 2003), firmly established. Higher mass would give higher 
C. We conclude that the value of constant C if likely to 
be between 1 and 2, though finding its more exact value 
would require thorough investigation of many X-ray binary 
sources. We plan to do this in the forthcoming paper. 

We showed that the two methods of black hole mass 
determination, with the scaling adopted as described, are 
statistically in agreement. The mean value of the mass ra- 
tios is r = 0.81 ± 0.15 (one sigma error), consistent with 1. 
Recent determination of the black hole mass in NGC 4151 
(4.14±0.73x 10 r M Q ) by Metzroth, Onken, Peterson (2006), 
based on reverberation, agrees even better with our values 
(based on excess variance) than the value given by Peterson 
et al. (2004) quoted in Table 1. This means that normal- 
ization factors used by the two methods were determined 
correctly, and their corresponding errors are perhaps even 
smaller than the conservatively estimated error. Therefore, 
black hole masses can be measured rather accurately, using 
any of the two methods if the requested data for an object 
is of appropriate quality. 

4.2 Mean geometrical factor 

Since more discussions in the literature usually concerned 
the squared geometrical factor / 2 , as introduced in equa- 
tion we convert for convenience the discussed results 
f 2 . The normalization obtained by Onken et al. (2004) corre- 
sponds to f 2 — 1.37T0.45. Our results from the X-ray excess 
variance method gave slightly lower values of the black hole 
masses so they effectively correspond to f 2 = 1.03 (median 
value) and 1.12 ± 0.20 (weighted mean). However, recently 
Collin et al. (2006) have obtained / 2 = 0.96 from the rever- 
beration method. 

Kaspi et al. (2000) use equation Q and adopt the value 
f 2 — 3/4 corresponding to a isotropic distribution of ran- 
domly oriented orbits (Netzer 1990). They also show that a 
convenient scaling exists between the continuum luminosity 
at A = 5100 A which is produced in accretion disc and the 
size of the BLR, -Rblr- This scaling was used later in many 
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papers (e.g. Woo & Urry 2002; Willot, McLure & Jarvis 
2003, Wu et al. 2004, Warner, Hamann & Dietrich 2004, 
Czerny, Rozariska & Kuraszkiewicz 2004) for estimation of 
the black hole mass without a need for time-consuming mon- 
itoring. Again, the factor f =3/4 was assumed. 

Results based on the excess variance support the con- 
clusion of Onken et al. (2004) that f 2 = 0.75 is certainly 
too low, and the appropriate value is by a factor ~ 1.7— 1.8 
higher than that. 

This will increase black hole masses derived with the 
use of Kaspi et al. (2000) formulae, and in consequence it 
will decrease the estimated value of the Eddington ratio. 
It will (at least partially) reduce the problem of the highly 
super-Eddington accretion rates in some quasars (Collin et 
al. 2002). 

Normalization higher by a factor of 2 than in the origi- 
nal formula of Kaspi et al. (2000) also leeds to estimates of 
the quasar radiation efficiency in better agreement with the 
timescales of the black hole growth (Wang et at. 2006). 



4.3 Broad Line Region Geometry 

Attempts of the direct determination of the BLR geometry 
from the reverberation studies of the line profiles in NGC 
5548 (Done & Krolik 1996) showed that none of the sim- 
ple models reproduced the velocity field of the clouds. It is 
possible that the Keplerian motion dominates but the disk 
dumpiness is responsible for the observed complexity of the 
picture (Shapovalova et al. 2004). 

Our weighted mean value of f 2 thus can be inverted 
into a mean value of the inclination angle, i, as done by Wu 
& Han (2001). If we neglect the random/wind component of 
the velocity (i.e. for £ = 0), mean value of f 2 corresponds 
to the inclination angle 28 ± 4°, and the median value of 
f 2 is equivalent to i ~ 23°. Both values are reasonable, 
taking into account that large viewing angles are obscured 
by the dusty /molecular torus and active nuclei with large 
inclinations are classified as type 2 objects. 



5 CONCLUSIONS 

• We determine the squared geometrical factor f 2 present 
in the reverberation method using black hole masses deter- 
mined from X-ray variance method as a reference. Our value 
of the geometrical factor, f 2 = 1.12, is consistent with the 
value 1.37 obtained by Onken et al. (2004) from the black 
hole mass - stellar dispersion relation, and applied by Pe- 
terson et al. (2004), and in agreement with f 2 = 0.96 of 
Colin et al. (2006). Such values are intermediate between 
the expected values for a spherical and a flat geometry of 
the BLR. 

• The values of the black hole masses given by the formu- 
lae of Kaspi et al. (2000) are most probably systematically 
underestimated by a factor of ~ 1.7-1.8. 

• complementary use of the reverberation and excess vari- 
ance methods is a promising tool to constrain the inclination 
of an AGN 
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